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Incompressibility, nonlinear deformation under stress and viscoelasticity are the fingerprint of soft 
tissue mechanical behavior. In order to model soft tissues appropriately, we must pursue to complete 
these requirements. In this work we revisited different soft tissue quasi-linear model possibilities in 
trying to achieve for this commitment. 



PACS numbers: J.1,J.2 
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I. INTRODUCTION 

Different mechanical properties of biological tissues de- 
pend on mineral content. Thus, roughly, it can be distin- 
guished two classes of biological tissues. Bones and tooth 
content mineral, they conform the group known as hard 
tissues. Whereas, skin, muscle, blood vessel and lungs 
conform the second group called of soft tissues. They do 
not content mineral, so they are much deformable than 
hard tissues. 

Modelling the mechanical behavior of soft tissues has 
much in common with those techniques used to model 
rubber elasticity. Thus, finite deformation theories useful 
for rubber elasticity are often used to describe soft tissue 
mechanical behavior. However, there are significant dif- 
ferences in the material structure of soft tissues and rub- 
ber elasticity. Moreover, in the way they respond under 
applied stress. Soft tissue material achieve initially large 
stretch under relatively low level of stress and subsequent 
stiffening at higher level extensions. This is shown in Fig- 
ure (EJ, where are compared the typical simple tension 
stress-stretch response of rubber (left-hand figure) with 
that of soft tissue (right-hand figure). For other side, col- 
lagen fibres distribution leads to pronounced anisotropy 
in soft tissue, in contrast with from typical isotropic rub- 

beriaaa. 

In biomechanics an specimen (animal) is regarded as 
an assemble of rigid bodies connected by joins and mus- 
cles. Under this framework, Newton's law balance of 
forces maybe could not be efficient to write equations 
of motion for each connectors, known their mechanical 
properties. Naturally, because the number of particles to 
be considered is excessively large, and the way to spec- 
ify the interaction forces between particles becomes very 
complex. Instead, it is much better to consider the body 
as a continuum, searching for a simplified way to specify 
interaction forces between particles by means of constitu- 
tive equations, and so thus specifying realistic properties 
of materials. A stress-strain relationship describes the 
mechanical property of a material and it is therefore a 
constitutive equation. 

Research on modelling soft tissue mechanical behavior 
has a growing demand for applications in surgical simu- 
lations, pursuing for in real time fast and precise calcu- 



lations of tissue deformations. In trying for these com- 
mitments, it has been introduced models accounting for 
the continuous nature of soft tissue. Within the limits of 
the employed constitutive model the finite element meth- 
ods allows for physically correct simulation of the tissue 
mechanics 0. 0. [t|. For other side, alternative discrete 
approaches based on spring-mass model have also been 
applied0,|3- It seems that methods based on the contin- 
uum approach are more appropriated to describe much 
better realistic soft tissue mechanical behavior, whereas 
spring-mass based models being economics are less accu- 
rate than the former. The main issue of our research look 
for an appropriate way to establish an accurate compar- 
ison between these methodologies. 

In this work we review some physic models applica- 
ble to modelling soft tissue mechanical behavior and we 
also discuss some perspective in the field. This work is 
organized as follow: In section |H] we present the contin- 
uum approach formalism, and we also consider particu- 
lar mathematical results for linear elastic solids and flu- 
ids. Important testing methods on physics necessaries to 
characterize soft tissue mechanical behavior, i.e, strain, 
creep and relaxation, are presented in section ITTT1 Next, 
in section II VI are reviewed more general soft tissue me- 
chanical properties. Section [V] is devoted to show basic 
mathematic models describing nonlinear elasticity and 
quasi-linear viscoelasticity, which are related to soft tis- 
sue mechanical behavior. In section IVT1 it is presented 
the basic vibrational spring-mass model, which is also ex- 
tended to consider viscoelastic mechanical behavior. In 
section IVTT1 we presented some discussions and perspec- 
tives. 



II. CONTINUUM APPROACH 

The equations of motion for a continuum were derived 
by Euler, using Newton's laws and supported on the fol- 
low axioms: i) The material particles form a one-to-one 
isomorphism with real numbers in a 3V Euclidean space. 
ii) The mass distribution is characterized by the density 
p (the mass per unit volume), it is defined as a piecewise 
continuous function over the volume of the continuum. 
in) Solely neighbors particle interactions are regarded. 
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FIG. 1: Typical simple tension response of (a) rubber and 
(b) soft tissue. Nominal stress t > plotted against stretch 
A > 1. 



A surface is conformed by particles, it could be consid- 
ered arbitrarily infinitesimal. In a continuum, oriented 
surfaces interaction can be expressed as a surface trac- 
tion ( Force per unit area), that can be computed from 
a well defined stress tensor (cr^ ). 

Let us xi,X2,xs be an inertial cartesian frame of ref- 
erence, and an infinitesimal volume by dx\dx<2.dx?,. Let 
cjjj, eij be the stress and strain tensor, respectively. Let 
us denote: Vi, Dvi/Dt and Xf, the velocity vector, the 
acceleration vector and the body force per unit volume, 
respectively. We now recall Newton's balance of forces, 
stating that the material rate of change of the linear mo- 
mentum of a body must be equal to the resultant of forces 
applied to the body. In follow, we write it in the form of 
Euler's equation, for more details see 0: 



This is valid for any continuum, whether it is a fluid or 
a solid. 

In particular, if the material is considered to be in- 
compressible, regarded p is a constant; we get from the 
Eqn.(E), 



dvj/dxj = . 



(4) 



From here, further development requires specification of 
the properties of the material in the form of constitutive 
equations, relating stress with strain or strain rate; or 
strain history. In order to have a complete theory, all 
these constitutive relationship must be incorporated into 
Eqn.(E|. 

Using our faculty perceptions we can verify: booth skin 
and muscle have a definite shape, they even can sustain a 
shear force, and maintains a quasi-static state (i.e. stress 
is a function of strain). Therefore, their mechanical be- 
havior is analogous to that one of a solid. Whereas con- 
sidering blood (at least in normal conditions) we note: it 
has no shape, it can not sustain a shear force (because it 
is in motion), and given its dynamic state; the stress is a 
function of the rate of strain. Thus it behaves mechani- 
cally like a fluid. 

Next, we derive particular equations for isotropic solid 
materials which follow Hook's elastic law and also for 
Newtonian fluids using Euler's formalism, incorporating 
specific constitutive relationships. 



A. Isotropic Hookean Solid 

These kind of materials obey the following stress-strain 
linear relationship: 



dp dpv-j 



Dvi da; 



Dt 



dxj 



+ Xi. 



(1) 



Using tensor notation, a repetition of an index in any 
given term means a summation over the range of that 
index, 



dxj 



dan da i2 da i3 



dxi 8x2 



dx:- 



In the spatial description, 
Dvi d v 



Dt 



dt 3 Ox, ' 
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(5) 



which alternative inverse relationship can be re-written 

as, 



1+y 
E 



1 1 &kk$ij i 



where a^ and are the stress and strain tensors, A is the 
Lame constant, G is the shear modulus, E is the Young's 
modulus and v is the Poisson's ratio. Sij is the Kronecker 
delta (unity Hi = j, otherwise is zero ). We see here that, 
the strain is being referred to a configuration of the body 
in which it is null provided that the stress is zero. 
Substituting Eqn.(EJ) into Eqn.(EJ) we arrived to, 



By taking into account the conservation of mass, ex- 
pressed by the equation of continuity: 



dp dpvj 
dt dxj 



Dp dvj 
Dt ^ dxj 



Eqn.(GJ becomes, 
Dvi dai 



P- 



Dt d Xj 



+ Xi 



0, 



1,2,3) 



(2) 



(3) 
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Dvi 
Dt 
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dxi 



dxj 



+ Xi 



(6) 



Lets now denote the vector displacement of a point in 
the body by Mj, which is measured with respect to an 
inertial cartesian frame of reference. If Ui(xi,X2,xs,t) is 
infinitesimal, we can define, 



1 / dui duj 

2 \ dxn dxj 



(7) 
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To the same order of approximation, 



duj 

dt ' 



Dvi 
Dt 



d 2 Uj 
dt 2 ' 



(8) 



On substituting all these equations into Eqn.( IHJ), it is 
obtained the well known Navier's equation: 



d 2 Uj 
' dt 2 



G 



d_Uj 
dxjdxj 



+ 



G d duj 



2u dxi dxj 



+ X t , (9) 



where the Poissons' ratio v 



In this way, we 



2(A + G)' 

arrived to the basic field equation of the linearized theory 
of elasticity. However, living bodies often take on finite 
deformations, obey so more complex nonlinear constitu- 
tive equations than the above presented. 



B. Newtonian Fluid 

These fluids obey the following stress-strain-rate rela- 
tionship, 



= -p Sij + [i 



dui duj 
dxj dxi 



\dx k 



&ij , (io) 



where p is the pressure, \x is the dynamic viscosity ( a 
constant parameter), and A is the bulk elasticity ( second 
coefficient of viscosity). If the fluid is incompressible, 
according to Eqn.(0J), this expression simplified to, 



-pSij + \x 



dui duj 
dxj dxi 



(11) 



By substituting it into Eqn.([2J), we obtain the Navier- 
Stokes equation, valid for Newtonian fluids, 



duj 
dt 



Uj dxj 



1 dp d 2 Ui 
pdxi dxjdxj 



+ X l7 



[12) 



where v (= \xj p) is the fluid kinematic viscosity. Appro- 
priate boundary conditions are necessary in order to solve 
the current equation. 

Air and water are Newtonian fluids. Non-Newtonian 
fluids do not obey constitutive equations ( Like 
blood, most body fluids are non-Newtonian. However 
at a shear strain rate above 100 sec -1 blood is almost 
Newtonian. 

In general, soft tissues present unique mechanical prop- 
erties. In the next section we present some basic test 
methods useful to describe quantitative soft tissue me- 
chanical properties. 



III. TEST METHODS 

The tensile test method consists to measure deforma- 
tion of an specimen while a force being applied to a sam- 
ple is gradually increased. During the stretching process 



are measured both, the stress and the strain. The former 
is obtained taking the ratio of the force and the material 
cross sectional area, whereas the latter is derived dividing 
the displacement and the original material length. Ana- 
lyzing the stress-strain curve provides useful information: 
as said the Young's module, the yield point and the ul- 
timate tensile strength. The Young's module is obtained 
taking the slope of the linear stress-strain curve. The 
yield strength corresponds to the value of the stress at 
the yield point, generically calculated by plotting Young's 
module at a specified percent of offset (0.2%). The ul- 
timate tensile strength corresponds to the highest value 
of stress on the stress-strain curve. Computing the area 
under the stress-strain curve we determine the fracture 
energy. 

Meanwhile, viscoelastic properties of soft tissues can be 
determined by creep and relaxation tests. Creep proper- 
ties are determined subjecting the specimen under pro- 
longed constant tension or compression, and the defor- 
mation is measured at specified times intervals. Then, a 
creep (deformation) versus time diagram is plotted, and 
the creep-rate is obtained by taking the slope of curve 
at any point. Provided failure occurs, the test is finished 
and the time for rupture is recorded. If the specimen 
does not fracture within the test period, creep recovery 
may be measured too. 

The stress-relaxation of a material is obtained taking 
the specimen deformed a given amount and recording 
the actual stress over prolonged period of exposure, and 
the stress-relaxation-rate is the slope of the curve at any 
point. 
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FIG. 2: Elastin-rich canine nuchal ligament, collagen-rich sole 
tendon, and intestinal smooth muscle tensile properties. The 
sole tendon shows to be harder than the other two kinds of 
tissues, smaller the deformation biggest is the strain. In con- 
trast, the intestinal smooth muscle shows to be more soft- 
ness, strain value of 0.5 is achieved with more facility than 
for nuchal ligament (elastin), it is also more viscoelastic than 
the others tissues, because it shows a wider hysteresis loop. 



IV. SOFT TISSUE BASIC PROPERTIES 
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Once the biologic specimen is subjected to test meth- 
ods as described above, their mechanical properties are 
obtained. In follow we describe more basic soft tissue 
mechanical properties. 



A. Inhomogeneous Structure 

Biological soft tissues are mainly composed of cells and 
intercellular substances. The later consisting of connec- 
tive tissues such as collagen, elastin and ground sub- 
stance (hydrophilic gel). These components present dif- 
ferent physical and chemical properties, and their con- 
tents differ from tissue to tissue and even from location 
to location within a tissue. Thus, mechanical properties 
depend both on tissue and on site. 

Collagen is a protein which gives mechanical integrity 
and strength to our bodies, it is presented in different 
structural forms in different tissues and organs. 

Elastin is also a protein, unlike to collagen has a much 
less integrated structure. It has much less strength and 
much more flexibility than collagen. 

Figure ( |2J) shows the stress-strain relationship for 
soft tissues rich in collagen, elastin and smooth muscle 
(cell) 10]. There we can see that, the elastin- rich nuchal 
ligament has much less strength and much more flexibil- 
ity than the collagen-rich sole tendon, whereas the in- 
testinal smooth muscle is much softer than the other soft 
tissues, which also shows a wider hysteresis loop; indicat- 
ing that it is viscoelastic. 
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FIG. 3: Tension-extension ratio relations of rabbit ventral 
skin. Different predetermined orthogonal directions are shown 
in the inset. 



Bovine articular cartilage 
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FIG. 4: Tensile properties of bovine articular cartilage for 
different orientations. Here, we can also note certain nonlinear 
strain response of the soft tissue under stress. 



B. Nonlinear Large Deformation and Anisotropy 

Biological soft tissues are mechanically nonlinear, this 
is enhanced by assembly of structural component into 
a tissue. The behavior of arterial walls at low tension 
is similar to that of elastin, while the behavior at high 
tension is the same as that of collagen. 

Collagen and elastin are long-chained polymers, so 
they are intrinsically anisotropic. In order to work ef- 
fectively, not only their fibers but also cells are oriented 
in tissues and organs. Skin for example has very different 
properties in two orthogonal directions. This is shown in 
Figure (|2Jl, where extensions under different applied val- 
ues of tension were measured. In this test method, Ten- 
sion consists to subject the specimen to pre-determined 
value of force, i.e. F X}V ] the subindex recall for the direc- 
tion, whereas the extension \ xy is measured, taking the 
ratio of material length under the actual tension and its 
original one (when the tension is null). Note here that 
the sample achieves it original length when the tension 
value is zero. 

Analogously, collagen fibers in the articular cartilage 
also shows anisotropic mechanical behavior 11] . This 
is presented in Figure ( HJ, where are shown different 
stress-strain curve responses for different orientations, 
i.e.0, 45, 90 degrees. 



C. Viscoelasticity 

Most rheologist would restrict using of the label non- 
Newtonian fluid to substances which, because they pos- 
sess long structural relaxation times, are non linear and 
tend to remember their past history. For almost all the 
biological soft tissues the curve stress-strain exhibit hys- 
teresis loop, as it was shown in Figure ( ^) for elastin, 
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FIG. 5: Relaxation behavior of rabbit patellar tendon 



collagen and soft muscle, meaning that they are viscoelas- 
tic. However, viscoelastic behavior is distinguished more 
clearly by relaxation test. For example, as it is shown 
in Figure ( [SJ , when the patellar tendon is elongated a 
given amount and the actual stress does not stay at a spe- 
cific level but decreases rather rapidly at first and then 
gradually^. 

Generically speaking, under different strain rates (test 
speed) biological soft tissues are mechanically not very 
sensitive. In addition, the area of the hysteresis loop 
does not depend on the strain rate. 



original cross sectional area specimen at zero stress). The 
most striking feature of this relationship is revealed in 
Figure ( EJ) when the slope of the stress-extension curve 
is plotted against T. We can fit the experimental data 
to a straight line, as first approximation, 



BT + C, 



(13) 



where A is the extension ratio, B and C are constants 
depending of the material. Prescribed T = 0, A = 1, 
from Eqn. ( I13|) we obtain, 



T 



§{exp[5(A 



1)] " 1} 



(14) 



Several other types of soft tissue such as the skin, the 
muscles, the ureter, the lung tissue, and articular carti- 
lage are found to follow similar relationship, Figure (|7J). 
So, it appears that the exponential type of stress-strain 
relationship is common to biological tissue. 

Realistically, in vivo tissues are not expose to such a 
pure tensile force conditions. Therefore, under more re- 
alistic mechanical conditions multi-axial test are often 
used to determine soft tissue mechanical behavior. So 
in several cases we need more equations. Naturally, for 
a three-dimensional organ we need a three-dimensional 
stress-strain law. 



B. Multi-axial Behavior 

No general constitutive equations has been identified 
for living tissues. Strain energy functions are used to 



D. Incompressibility 

Most biological tissues hardly change their vol- 
ume even if load is applied, so they are almost 
incompressible [l3j. They have a water content of more 
than 70 %. However, articular cartilage is an exception, 
because is a micro-porous tissue, therefore water can en- 
ter an leave pores depending upon load 14]. 

The incompressibility assumption is a very impor- 
tant ingredient in constitutive laws formulation, impos- 
ing that adding overall principal strains is always zero. 

V. MATHEMATICAL DESCRIPTION OF 
MECHANICAL BEHAVIOR 

A. Tensile Behavior 

Biological soft tissue can support large deformations 
which are nonlinear and extensive. In reducing the ex- 
perimental data to a stress-strain relationship, it is used 
the nominal stress T (getting dividing the force by the 
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FIG. 6: Slope of stress-extension ratio curve versus stress of 
human dura matter and bovine pericardium 
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direction is given by 
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FIG. 7: Stress-extension curves of human dura matter and 
bovine pericardium . 



describe the mechanical tissue properties. The main idea 
is that if a pseudo-elastic strain energy function exist, 
then the stress-strain relationship can be obtained by a 
differentiation. Let W be the strain energy per unit mass 
of a tissue, and po be the density in the zero-stress state. 
Thus po W is the strain energy per unit volume of a tissue 
in the zero-stress state, this is called the strain-energy 
density function. The strain energy W can be expressed 
in terms of the components Eij of the Green strain. In 
a uni-axial tensile test, the Green strain E in the tensile 



E — E\-\ 



(15) 



The components Sij of the second Piola-Kirchhoff stress 
can be obtained by derivatives of the strain energy den- 
sity function, 



S; 



d Po W 
~dE~ 



(16) 



In a simple uniaxial tensile test, the only non-zero com- 
ponent of stress, denoted by S, is 



S = S U 



F 
~\A( 



(17) 



where Aq and F are the initial cross-sectional area sample 
and the load applied to the specimen, respectively. 

Equation ( ITrjj) is a constitutive law or stress-strain 
relationship. Thus we need to know the details of the 
strain-energy functions. For the strain energy function, 
studies on non isotropic tissues as skin, muscles, liga- 
ments, etc., have shown that the exponential form ap- 
plies well. A general strain-energy function for biological 
soft tissues is given by Fung 9] in the form 



W 



i&ijkiEijEk 



x exp {^iijEij + jijkiEijEki + . . .) , 



E E \ 



(18) 



where acjki, Po, Pmnpq, Jij,- ■ ■ are constants. 

Regarded the complex form of Eqn. ( I18|) high order 
terms can be dismissed as for skin, considered as being 
in a bi-axial state, and its strain-energy function is, 



po W = § [pt x El x + a 2 E$ 2 + a 3 E^ 2 + a 3 E% x + 2a 4 E 11 E 22 ) 

\cexp ( ai E^ + a 2 E% 2 + a 3 E\ 2 + a 3 E^ x + 2a A E xl E 22 + lx Ef x + l2 E% 2 + <y 4 E? x E 22 + 75^11^22) » 

I 



(19) 



where ai, . . . , 04, a±, . . . CZ4, 71, . . . , 75 and c are con- 
stants. 

Following, using Eqn. ( TTB|). we obtain the stresses, 



^- = a 1 E 11 + a 4 E 22 +cA 1 X, 



Sn 

S22 =^ = a 4 E 11 +a 2 E 22 +cA 2 X, (20) 



'12 - -m^ 



a 3 Ei 2 + ca 3 E\ 2 X , 



where A±, A 2 and X are functions of the strain compo- 



nents. 

It can be verified from Figure ( |3J) that Eqn. ( I19|) is 
suitable for bi-axial mechanical behavior of rabbit skin. 
Additional terms can be disregarded taking in account 
physical considerations, if we are concerned mainly with 
higher stresses and strains the first group of terms in Eqn. 
( I19|) can be omitted, in addition for practical purposes 
third order terms in the exponential function can also be 
omitted, and we have simply, 



Po W = \cexp (aiE xl + a 2 E 22 + a 3 E X2 + a 3 E 21 + 2a 4 E 11 E 22 ) . 



(21) 
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C. Viscoelastic Behavior 

Hysteresis, relaxation, and creep are common features 
of viscoelastic behavior for many materials. The toy 
models here presented are basic mechanical models useful 
describing viscoelastic behavior. These are composed by 
a combination of elastic material with elastic constant fi 
(i.e. spring) and a dashpot containing a fluid with coef- 
ficient of viscosity r\. A force F is suppose to produce on 
the spring an uniaxial linear deformation, denoted by the 
displacement u, and proportional to \x. Whereas, at any 
instant, for the dashpot, it produces a velocity propor- 
tional to r\. Different spring-dashpot arrangements are 
shown in Figure ([HJ). 

It is remarkable that these classes of discrete mod- 
els are also useful in the continuum, where relaxation 
functions are borrowed from these theories. For exam- 
ple modelling viscoelastic tissue behavior is assumed that 
the stresses in the material, which themselves may result 
from a nonlinear stress-strain relation are linearly super- 
posed with respect to time, recalling the quasi-linear vis- 
coelastic formulation, which we present next in subsec- 
tion (En). 



1. Maxwell Model 



sudden application of a constant unitary force at t = 0, 

c(t)= Q + V)l(t). (23) 

Interchanging the role of F and u, we obtain the relax- 
ation function g(t), as a response to a unity elongation, 
i.e. u(t) = l(t). The relaxation function is the force that 
must be applied to the specimen in order to produce an 
elongation, that changes at t = from zero to unity, 
remaining unity thereafter. 

g(t)=tiexp(-Ht/v)Ht), (24) 
with, l(t) is the step function, 

( 1 ; t>0, 
l(t) = I 1/2 ; t = 0, (25) 
( ; t<0. 

In this way, for the Maxwell model we found that a sud- 
den application of a load induce an immediate deflection 
of the elastic spring, followed by a creep of the dashpot, as 
in Eqn. ( I23J) . On the other hand, a sudden deformation 
produces an immediate reaction by the spring followed 
by stress relaxation according to exponential law, which 
relaxation time is given by the factor rj/^t, as in Eqn. 

(E§. 



In this model an spring and a dashpot are arranged in 
serie. Thus, the material behaves like an elastic at short 
time and it is viscous at long times. 

u = F/fi + F/ V ; «(0) = ™ • (22) 

Equation ( I22J) is solved, prescribed the force F(t) is a 
unit step function, and the result is called the creep func- 
tion c(t). It represents the elongation produced by a 
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FIG. 8: Spring-dashpot arranged as to describe different basic 
viscoelastic models. 



2. Voigt Model 

Let us now regard the spring and the dashpot in par- 
allel, 

F = nu + rju ; «(0) =0. (26) 

Solving Eqn. ( I2rjj) we obtain the creep and the stress 
relaxations functions, respectively, 

c(t) = Ml-exp (-/it/77)} l(t), (97 , 

g (t) = nS(t) + iii(t). [zn 

Where, 5(t) is the Dirac delta function. Regarding Eqn. 
( I27|) we found that, a sudden application of a unity 
force will produce no immediate spring deflections; which 
solely built up gradually, because the dashpot arranged 
in parallel with the spring will not move instantaneously. 
The displacement of the dashpot will relax exponentially 
with a relaxation time rj/fj,. 

3. Kelvin Model 

This model is also known as the Standard Linear 
model. While models due to Maxwell and Voigt suppose 
fluids behave like elastic to some extent, the standard 
model goes further taking in account the dissipation rate 
of energy in various materials subjected to cyclic loading, 

F + r e F = E R (u + r a u) , r e F(0) = E R r a u(0) .(28) 
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Here r e is the time of relaxation of the load under the 
condition of constant deflection, whereas r a is the time 
of relaxation of deflection under the condition of constant 
load, and 



c(t) 

g(t) 



Evil 



1 - £ 
1 - ^ 



-t/7 
-t/l 



l(t). 



(29) 



In Eqn. ( I29j) as t — ► oo, the dashpot is completely 
relaxed and the load- deflection relaxation becomes that 
of the springs, characterized by the constant Er, which 
is called the relaxed elastic modulus. 



4- Boltzmann Model 

This is a most general formulation under the assump- 
tion of linearity between cause and effect. In the 1 — T> 
case, analogously to the case above, we consider a bar 
subjected to a force F(t) and an elongation u(t). The 
elongation is caused by the total history of the loading 
up to the time t. Lets now consider an small time inter- 
val dr; at time r, the increment of loading is (jj^) dr . It 
contributes an element du(i) to the elongation at time t, 
with a proportionality constant c, depending on the time 
interval t — r , then, 



f* / v dF 
u(t) = / c(t — t) —— dr . 
Jo dr 

Analogously, for the force, 

\ f / \ du(r) , 
F(t) = J g(t-T)-±±dT. 



(30) 



(31) 



Here the functions c(t—r) and g(t—r) are the creep func- 
tion and the relaxation function, respectively. More gen- 
erally, we can write the relaxation function in the form, 



N 



g(t) = ^2a n e~ 



(32) 



n=0 



where v n = l/r n is called characteristic frequency and 
a n is an amplitude associated to each frequency v n , it is 
called a spectrum of the relaxation function. A gener- 
alization to a continuum spectrum is giving in the next 
section. 



D. Quasi-Linear Viscoelasticity of Biological 
Tissues 



The hysteresis H is due to viscoelasticity, it is defined 
as the ratio of the area of the hysteresis loop divided by 
the area under the loading curve. For biological tissues, 
H is seen to be variable, but its variation with strain 
rate is not large, see Section f llV CJ) . Hence, typical soft 
tissue behavior shows non-linearity of the stress-strain 



relationship and insensitivity of the material to strain 
rate. 

Models due to Maxwell and Voigt provide a logarith- 
mic type functional relationship between hysteresis and 
frequency (v, the inverse of the characteristic time), i.e., 
H{u) oc ln(z/). It is a decreasing curve for the former, 
whereas it is a increasing curve for the later. On the other 
hand, the Standard model presents a bell-shaped curve 
of hysteresis vs the logarithm of the frequency. Summa- 
rizing, none of these models present a typical flat hys- 
teresis vs frequency curve of living tissues. However, this 
fact can be corrected by introducing nonlinear springs. 
Hence, this now suitable model for soft tissues, is com- 
posed of a long series of Kelvin bodies whose character- 
istic times (l/v) span a broad range. 

Lets consider an elastic stress tensor T^ e \ which is a 
function of the strain tensor E. If the material is in the 
zero-stress state until the time t = 0, then it is suddenly 
strained to E and maintained constant at that value, thus 
the stress developed will be a function of time as well as 
of E. Hence, the history of the stress development may 
be written as, 

Gijmnit) ^mn(^) j Gij mn (0) = 1 j (^3) 

in which Gij mn (t) is called the relaxation function, a nor- 
malized function of time. Now, we suppose that the stress 
response to an infinitesimal changes in a component of 
strain, SEij, superposed on a specimen in a state of strain 
E at an instant of time r, is, for t > r , 



(t-r) 



dT£l[E(r)] 
dEk* 



dE l3 {r) 



(34) 



In addition, assuming that the stress at time t is the sum 
of the contributions of all the past changes, each governed 
by the same reduced relaxation function, 



Tij{t) 



dT^l[E(r)] dE kl (r) 



G 



ijmn 



(t-r) 



dE kl 



dr 



dr .(35) 



Although T^(E) may be a nonlinear function of strain 
the relaxation process is linear. Hence, this is a quasi- 
linear viscoelastic theory. 

In particular, for the one-dimensional Kelvin model, 

l + S(q)exp(-t/q) 



G(t) 



(36) 



1 + S(q) 

with q and S(q) are constant parameters. The latter is 
called the relaxation spectrum and \ jq is a frequency. 
For an infinite number of Kelvin model in series, we can 
get the follow reduced relaxation function in the form, 

l + / °° S(q)exp(-t/q)dq 



G(t) 



l + J OO S(q)dq 



(37) 



It has been shown that a specific spectrum, with con- 
stants c, qi, q2, the follow function fits the data for the 
skin, arteries, ureter and teniae coli. 



S{q) 



c/q ; qi< q< q2, 
; q < q 1 ; q > q 2 . 



(38) 
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VI. VIBRATIONAL SYSTEMS 

In this section we suppose that a physical body is con- 
formed by an ensemble of connected springs. Initially is 
being considered both spring elasticity and spring dump- 
ing responses are linear. 

Hence, let us consider an elastic body obeying Hooke's 
law. Let it be rigidly supported in a fixed space and sub- 
jected to a set of forces Qi acting at points i £ (1, . . . , n). 
Lets set up a generalized displacement at point i, in 
the direction of the force Qi, be linearly proportional to 
the forces Qi, . . . , Q n , and vice versa, for i £ (1, . . . , n): 



q> 
Qi 



CnQi +Q2Q2 
Knqx +K i2 q2 



+ ... 
+ ... 



Ci n Q n , 
Kin Qn ■ 



(39) 



The constants of proportionality Cij , Kij are indepen- 
dents of the forces and displacements. Cij are the flex- 
ibility influence coefficients and are the stiffness in- 
fluence coefficients, respectively. The physical meaning 
of Kij is the force required to act at the point % due to 
a unit displacement at the point j, while another points 
(others than j) are held fixed. Analogously, Cij is the 
deflection at i due to a unit force acting at j. The con- 
straining equations ( I39|) imply that exist a unique un- 
stressed state to which the body return whenever all the 
external forces are removed, so the superposition princi- 
ple applies. Moreover, it implies that the total work done 
by a set of forces does not depend of the order in which 
the forces are applied. This work done by each force is 
\QiQi and the total work done by the system of forces, 
is stored as strain energy in the elastic body , and it is 
given by, 



u = \ ElLi QiQi = \ £"=1 Ej=i KijQi9j , 

= 2 Ei=i Ej=i Qi Qj ■ 



(40) 



The matrices Kij and Cij are symmetric. In others 
words, the displacement at a point i due to a unit force 
acting at another point j is equal to the displacement at 
j due to a unit force acting at i.e., they are positives 
in the same direction at each point. The stability of the 
system is guarantied if all the principal minors, including 
the determinant of the full Kij {Cij) matrices are posi- 
tives. 

Now we can write the equation of motions of a set o 
masses embedded in an elastic body. The D'Alambert's 
principle establishes that the particles can be considered 
to be in a state of equilibrium if the inertial forces are 
applied in reversed direction on the particles. Thus, if the 
body is attached to a stationary support and we consider 
dumpers (lets said, dumping forces), then the dynamic 
equations written in vectorial form are 

n n 

mi'qi + ^2 K H Qj + z2^ijQj = F i ! (» = 1» • • • n ) ■ ( 41 ) 
j=i i=i 

Here qi is the displacement of the particle and Ff is the 
external force acting on it. 



However, damping force could be much more complex 
than this, for example it can be aerodynamics in origin, 
nonlinearly viscoelastic (as for tissues) and not necessar- 
ily stabilizing. For an immediate generalization of the 
theory, let recall for the concept of kinetic energy /C, it 
is an homogeneous quadratic function of the q's, and its 
rate can be computed from, 



x ^ / d dlC dJC\ . 
= ^\dtdq-~dq-) qt 



(42) 



For a biological system, other forms of energy rather 
than the kinetic energy K. can also be involved, such as 
the gravitational potential Q, the internal energy IA, and 
the chemical energy C. According to the first law of 
thermodynamics, the energy of a system solely can be 
changed by absorption of heat TL and by doing work on 
the system. The rate of change of the total energy must 
be equal to the sum of the rates of heat input H and 
work done on the system VV, so, 



ic + g + u + c = H + w. 



(43) 



Generally, just a certain part of the internal energy func- 
tion U arise from the strain energy, it is a quadratic func- 
tion of (generalized) coordinates qi. It can also depend on 
the temperature. If the temperature remains constant, 
then the rate at which work is done by the internal energy 
and also by the forces acting on a system are, respectively, 



then, substituting them into Eqn. ( H3|) . 



(44) 



E 

i=l 



d_d1C 
dt dq\ 



dlC dU 

dqi dqi 



Qi)qi=n-G-C. (45) 



Now lets consider the special case when H — Q — C 
f) . Because the displacements are independent variables, 
that can assume arbitrary variables, we impose for any 
qi / 0, in Eqn. ( H5Jl , 



d dlC dlC + dU _ 
dtdqi dqi dqi 



(46) 



If the (generalized) forces are partial derivatives of a func- 
tion, i.e., Qi = — jPp; where W is the potential energy. 
For the spring body system the internal elastic energy is 
given by Equation f lUIJ) . Thought, W and U only depend 
on qi, the system is conservative. 

The aerodynamic or hydrodynamics forces acting on 
an animal moving in a fluid are not conservative. Hence, 
for problems involving these forces is much better to use 
Equation (EBJ), thought Qi is an external force. 

Following, we apply the present theory to materials 
showing viscoelastic behavior (like tissues), in the pres- 
ence of fluid dynamic forces. 



10 



A. Systems with Dumping and Fluid Dynamics 
Loads 

Here we recall the most general constitutive equation 
of a linear viscoelastic solid reviewed in Section IIV CI 
Applying Boltzmann's formulation we may write the in- 
ternal force acting on a particle i, located at point i, by 
a displacement qj located at point j as, 

F tj (t) = J* G^t-^^^dr, (47) 

where Gij(t) are the relaxation functions. For solid with 
fading memory, it is given by equation 

Gij{t) = Kij + 0ije- tVl + ... + Hje- tv - , (48) 

where V\,V2, ... , v n are the relaxation frequencies , 
K^, Pijjjij are spectral constants. 

For systems subjected to fluid dynamics forces (i.e., in 
swimming, in wind or with internal flow), the external 
force at the point i due to motion at point j may be 
written as, 

F^(t) = f^A tJ (t-r)^dr, (49) 

where Aij(t) are aerodynamic influence functions, which 
is un-symmetric. 

Considering both the internal forces and a fluid dy- 
namic forces given by Ean.f 147)1 and Eqn. ( 149)1 . Eqn. 
( EE)) can be rewritten as, 

m, d -^ + T.UJLG„(t-r) d -^dr 

where F± is the external force other than fluid dynamic 
forces. 

VII. CONCLUSIONS AND PERSPECTIVES 

More general mechanical properties of a material can 
be specified by mean of constitutive equations. Remark- 



able, the continuum approach naturally adds constitu- 
tive equations into more generalized theories. However, 
it seems that soft tissue models based on these more gen- 
eral theories are numerically expensive, hence that linear 
elasticity approaches have been preferentially considered. 

On the other hand, spring-mass models present an- 
other alternative to model soft tissue mechanical behav- 
ior. Because of physical model simplicity and easier im- 
plementation, it has been preferentially used on applica- 
tions for surgical simulations. The method is economic 
and fast, so it is a good candidate to operate in real 
time. However, until it is of our knowledge, its numeri- 
cal implementation has been done no far away from the 
limits of linearity, not only regarding elasticity but even 
viscoelasticty. 

In this work we have review some possible alternative 
models considering nonlinear elasticity and quasi-linear 
viscoelasticity on spring-mass models. The former can 
be naturally filling up by adding to the model exponen- 
tial type deformations, as in section IV Al Viscoelastic 
mechanical behavior can be added to the model using 
quasi linear models based on Boltzmann type approach, 
as shown in section IVI Al Summarizing, these are the 
formulations we are going to explore in order to better 
model soft tissues. The resulting model will be validate 
by comparison with accurate models based on the contin- 
uum approach. We shall stress that because of simplicity 
and economy, this more elaborated spring-mass models 
can still handle to perform real time calculations. 

Recently, Holzapfeld [T(| and Ogden [l7| have inves- 
tigated numerical method based on the Continuum ap- 
proach in order to regard nonlinear elasticity and Vis- 
coelastic response on arterial walls and soft tissues. In 
addition to research for more sophisticated spring-mass 
models, these numerical methods are also going to pave 
the way of our research. 
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